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ABSTRACT 

Recent  advances  in  spiral  bevel  gear  geometry  and  finite 
element  technology  make  it  practical  to  conduct  a  structural 
anatysis  and  anatytically  roll  the  gearset  through  mesh.  With  the 
adv^  of  user  specific  programming  linked  to  3D  solid  modelers 
and  mesh  g^serators,  model  generation  has  become  greatly 
automated.  Contact  algoriduns  available  in  general  purpose  finite 
dement  codes  eliminate  the  need  for  the  use  and  alignment  of  gap 
elements.  Once  the  gearset  is  placed  in  mesh,  user  subroutines 
attadied  to  the  FE  code  easily  roll  the  gearset  through  mesh.  The 
method  is  described  in  ddail^  Preliminary  results  for  a  gearset 
segment  showing  the  progression  of  the  contact  lineload  is  given 
as  the  gears  roll  through  mesh. 

INTRODUCTION 

One  of  the  most  conq)lex  and  highly  loaded  components  in  a 
helicopter  is  the  main  rotor  transmission,  (Cormier  1979).  One 
of  the  critical  compon^its  of  the  transmission  is  the  spiral  bevel 
gear  mesh. 

The  designer  is  faced  with  tradeoffs  between  weight  and 
reliability.  Weight  reduction  results  in  increased  flexibility  and 
deflection  of  the  teeth  and  rim.  Deflection  of  the  gear  creates  a 
departure  from  the  expected  contact  found  by  rolling  the  gears. 
These  deflections  cause  the  contact  zone  to  shift.  The  shift  in 
contact  can  dramatically  increase  contact  and  tooth  bending 
stresses  with  a  corresponding  reduction  of  component  life. 

The  contact  elL^^se  and  localized  stress  distribution  have  been 
historically  analyzed  with  Hertzian  theory  developed  for  ball  and 
roller  bearing  elements  (Coy  et  al.  1985).  This  theory  considers 
the  contact  stresses  for  two  general  surfaces  with  direct  bearing 


load.  Gear  teeth  are  not  loaded  directly.  The  load  is  transmitted 
through  flexible  teeth  coimected  to  flexible  rims.  The  tooth  and 
rim  flexing  of  the  gear  and  pinion  will  greatly  effect  the  location 
and  shape  of  the  contact  zone,  load  sharing  between  teeth  and  the 
contact  stress  distribution.  These  factors  will  in  turn  affect  fatigue 
life,  noise  and  vibration,  thermal  degradation  and  hydrodynamic 
action  of  the  lubrication.  These  issues  are  also  important  design 
parameters  for  spiral  bevel  gearing  in  ground  vehicles. 

Considerable  work  has  been  completed  to  automate  3D 
analysis  of  spiral  bevel  gears.  Fundamental  to  this  work  is  an 
understanding  of  the  gear  surface  geometry.  Litvin  (1989)  and 
Litvin  et  al  (1989)  presented  the  theory  of  spiral  bevel  gear 
generation  and  design.  Handschuh  and  Litvin  (1991)  extended  this 
anatysis  to  the  math^natical  description  of  the  gear  tooth  surfaces. 

Bibel,  et  al  (1993, 1994)  began  die  automation  of  the  analysis 
with  model  development  software  that  creates  input  for  the  3D 
solid  modeler  (PATRAN)  and  the  finite  elem^  code  NASTRAN. 
The  finite  element  modeling  done  with  NASTRAN  utilized  gap 
elements.  The  analysis  was  greatty  sinqilified  by  using  the  contact 
algorithm  in  the  finite  element  code  MARC,  Bibel,  et  al  (1994). 
A  user  subroutine  compiled  with  the  MARC  input  file  rotates  the 
gearset  through  mesh. 

This  work  is  being  extended  into  a  series  of  programs  that 
receive  machine  tool  settings  as  input  and  create  3D  contact 
analysis  spiral  bevel  gears  rolling  through  mesh. 

SURFACE  GEOMETRY 

The  most  difficult  task  associated  with  3D  modeling  of  spiral 
bevel  gears  in  mesh  is  an  accurate  description  of  the  gear  surface 
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geometiy.  The  surface  of  a  generated  gear  is  an  envelope  to  the 
family  of  surfaces  of  the  head  cutter.  In  other  words,  the  points 
on  the  geaerated  tooth  surface  are  points  of  tangency  to  the  cutter 
sur&ce  during  manufacture. 

The  conditions  necessary  for  envelope  existence  are  given 
kinematically  by  the  equation  of  meshing.  The  equation  is  stated 
as  follows:  the  normal  of  the  generating  surface  must  be 
perpendicular  to  the  relative  velocity  between  the  cutter  and  the 
gear-tooth  surface.  This  is  given  mathematically  as: 

n-V  =  0 

Where  n  is  the  normal  vector  and  V  is  the  relative  velocity. 

The  equation  of  meshing  for  straight-sided  cutters  with  a 
constant  ratio  of  roll  between  the  cutter  and  worlqpiece  is  given  in, 
Litvin  et  al,  (1989  and  1991)  as: 

(u  -  r  coti(r  cosi|r)  cosy  sinx  -h 

S[(nicw  “  siuy)  cosT|r  sin0  t  cosy  sim|r  sin(q  -  <|>  J] 

X  E„(cosY  sinijf  +  siny  cosTjr  cost)  - 

siny  cosijr  sinx  =  0  Eq.  1 


Where:  X  = 
Y  = 

u,  0 

♦  = 
r  = 

= 


q  = 


(0  T  q  ±  ) 

the  root  angle  of  the  component  being 
manufactured. 

locate  a  point  on  the  cutting  head  as  given  in 

(Handschuh  et  al  1991). 

the  cradle  ori^itation 

the  blade  angle 

the  radius  of  the  cutter 

die  vector  sum,  a  machine  tool  setting  described 
in  (Handschuh  et  al  1991). 

M^hine  offset,  a  madune  tool  setting  described 
in  (Handschuh  et  al  1991). 
cradle  to  cutter  distance 
cradle  angle 

ratio  of  angular  velocity  of  the  cradle  to  the 
workpiece. 


All  of  the  above  terms,  except  u,  0  and  are  known 
machine  tool  settings  for  a  given  gearset  design.  These  three 
terms  can  be  mimericalty  solved  for  if  two  additional  equations  are 
identified. 

The  points  on  the  gear  surfaces  fouixl  by  solving  the  equation 
of  meshing  must  also  satisfy  the  axial  and  radial  location  of  a 
projection  into  the  XZ  plane  that  identifies  the  generated 
workpiece.  This  is  satisfied  by  the  following: 

Z;^-Z'  =  0  Eq.2 

r'-(X,*  + Y.T*  =  0  Eq.3 


Where:  and  are  coordinates  of  a  point  on  the 

worlqpiece  in  the  workpiece  coordinate  system,  and  Z*  and  r'  are 
the  a:^  and  radial  projections  into  the  XZ  plane  of  the  woriq)iece 
as  shown  in  Fig.  1. 

Equations  1, 2  and  3  can  be  solved  with  a  numerical  iterative  v 
procedure.  Equation  1  is  given  in  terms  of  cutting  head  coordi-  ^ 
nates  u,  0  and  (|)^ ,  whereas  equations  2  and  3  are  given  in  terms 
of  workpiece  coordinates  Xi^Y^and  Z^*  • 

A  point  on  the  head  cutter  (written  in  terms  of  u,  0  and 

must  be  transformed  into  a  point  on  the  workpiece  (written  in 
terms  of  X^Y^and  Zy^)  with  the  following  series  of  homogeneous 
coordinate  transformations. 

r,  =  IMJDVy  [M^JpvUIMJr,  Eq.  4 

Where:  [M  J  [M  J  [MpJ  [M  J  and  [M  J  are  a  series  of 
matrix  coordinate  transform^ions  given  in,  Litvin  (1989),  Litvin 
et  al  (1989),  and  Handschuh  et  al  (1991). 

These  coordinate  syst^ns  are  described  as  follows  and  shown 
in  Fig.  2  for  a  right4iand  gear  surface.  The  head-cutter  coordinate 
system  S^,  is  rigidly  connected  to  the  cradle  coordinate  systems 
S,.  rotates  about  fixed  coordinate  system  S„  attached  to  the 
fnflnbinft  fiiBine.  Coordinate  syston  Sp  orients  the  pitch  apex  of  the 
gear  being  manufactured.  The  common  origin  for  coordmate 
^sterns  Sp  and  locates  the  apex  of  the  gear  under  consideration 
with  respect  to  system  Sj^.  The  final  transformation  is  from 
coordinate  system  to  which  is  fixed  to  the  component  being 
manufactured. 

A  computer  program  was  written  using  the  above  solution 
technique.  This  program  receives  machine  tool  settings  and  basic 
gear  geometry  as  input  and  creates  an  output  file  containing  the 
coordinates  of  surface  points  for  the  gear  and  pinion.  The  output 
from  this  program  is  read  by  a  second  computer  program.  This 
program  uses  tibe  surfu^e  points  as  input  and  creates  an  output  file 
that  can  be  conpiled  and  read  into  the  3D  solid  rxKxleler  PATRAN 
(1991).  PATRAN  is  used  to  create  a  NASTRAN  or  MARC  input 
deck  of  any  mesh  density. 

A  gearset  model  created  in  this  manner  is  shown  in  Fig.  3. 
The  design  data  is  given  in  Table  I.  The  description  of  these  two 
programs  and  coding  is  given  in,  Bibel  et  al.  (1993  and  1994). 


GAP  ELEMENTS 

Previously  3D  contact  analysis  was  accomplished  using  gap 
elements.  This  presented  numerous  problems  related  to  gap 
element  orientation.  After  the  contact  point  was  located,  the 
surfece  normal  at  the  contact  point  was  calculated.  Lines  parallel  a 
to  the  contact  point  surface  normal  were  constructed  from  nodes 
on  the  pinion. 

The  intersection  of  these  lines  with  the  gear  surface  was  then  v 
calculated.  These  intersection  points  on  the  gear’s  surface  were 
then  incorporated  into  new,  distorted  elements  on  the  gear.  For 
every  new  gear  surface  node  that  was  added,  an  existing  gear 
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surface  node  had  to  be  deleted  and  the  element  nodal  connectivity 
had  to  be  redefined.  For  coarse  models,  it  worked  reasonably 
well  to  select  the  new  element  connectivity  by  visual  inspection. 
Larger  models  would  require  calculation  of  various  permutations 
of  nodal  connectivity  and  distances  to  determine  the  most 
appropriate  nodal  substitutions  arid  new  element  connectivity.  The 
resultiog  remiq^ped  elements  on  the  gear  surface  are  shown  in  Fig. 
4. 

This  procedure  worked  reasonably  well  for  coarse  models. 
However,  to  roll  the  gearset  through  mesh,  the  entire  process  of: 
contact  point  and  contact  point  surface  normal  calculation, 
iitfersection  of  lines  from  pinion  nodal  points  with  a  gear  surface, 
definition  of  new  element  connectivity  on  gear  surface,  and 
calculation  of  gaps  had  to  be  repeated  for  every  incremental 
rotation.  Typical  outyut  from  a  gap  element  model  is  shown  in 
Fig.  S.  Although  the  model  was  relatively  coarse,  the  results 
compared  favorably  with  Hertzian  estimates. 

CONTACT  ALGORITHM 

Greatly  superior  to  the  method  of  gap  elements  is  the 
automated  contact  algorithm  described  in  MARC  K^S.l  (1981). 
Contact  between  deformable  bodies  is  handled  by  imposing  non> 
penetrating  constraints  as  shown  in  Fig.6.  Automatic  in  this 
context  means  that  user  interaction  is  not  required  in  treating 
2iBilti>body  contact.  The  program  has  automated  the  inqxjsition  of 
iKn^paietration  cmistiaints,  evoi  while  rolling  through  mesh.  The 
user  no  longer  has  to  worry  about  the  location,  orientation,  arrd 
open/close  status  checks  of  “gq>  elements.**  Typical  output  for 
contact  analysis  from  the  contact  algorithm  is  shown  in  Fig.  7. 
The  results  were  consistent  with  previous  modeling  done  with  gap 
elements  and  with  Hertzian  estimates. 

AUTOMATED  MESHING  WITH  USER  SUBROUTINE 

The  model  development  software  nrathematically  creates  the 
pinion  and  gear  surface  coordinates  along  the  same  axis  of 
rotation.  The  gear  and  pinion  must  be  properly  located  for 
meshing  action.  First  the  pinion  is  rotated,  a  small  angle  about 
its  axis  of  rotation,  to  accommodate  tooth  mesh  alignment. 

The  jnnion  is  i^aced  in  mesh  by  a  rotation  6  corresponding  to 
the  shaft  angle.  The  gearset  is  now  properly  positioned  for  3D 
contact  analysis  with  MARC. 

After  corxq>letion  of  an  increment  of  the  finite  element  contact 
analysis  with  MARC,  the  gear  and  pinion  are  rotated  by  a  user 
subroutine  to  simulate  the  gears  rolling  through  mesh.  For 
the  gear  is  rotated  a  small  increment  of  “6”  degrees.  The 
pinion  would  then  rotate  “6  •  (Number  of  Gear  Teetii)/(Number 
of  Pinion  Teeth)”  in  the  opposite  sense.  The  finite  element 
contact  analysis  is  repeated  for  a  new  orientation  of  the  gearset. 
The  process  is  repeated  by  merely  changing  the  two  numbers 
corresponding  to  the  rotations  of  the  gearset  in  the  user  subroutine. 


The  user  subroutine  is  written  in  FORTRAN  and  compiled 
with  the  MARC  input  file. 

PRELIMINARY  RESULTS 

A  two-tooth  pinion  segment  and  a  two-tooth  gear  segment 
were  rolled  through  mesh  with  increments  of  about  one  degree  of 
pinion  rotation.  (There  are  19  teeth  on  the  pinion.) 

A  series  ofplots  are  shown  in  Fig.  8.  Shown  is  a  plot  of  line 
loading  of  the  forces  on  the  pinion  as  it  rolls  through  mesh  in  a 
counter  clockwise  direction.  The  line  load  is  seen  to  progress 
from  the  heel,  of  the  back  pinion  tooth,  to  the  toe. 

The  loading  switches  from  the  back  pinion  tooth  to  the  front 
afi^  a  brief  period  of  load  sharing.  Also,  shown  is  an  example  of 
line  loading  on  the  gear  for  a  single  configuration  in  Fig.  9. 

More  woxk  remains  to  be  done  concerning  mesh  refinement, 
enhanced  accuracy  of  the  modeling  of  the  fillet  region,  better 
boundary  conditions,  etc.  The  work  at  this  point  is  considered 
analogous  to  initial  attempts  at  2D  FE  modeling  of  spur  gears. 
The  models  went  through  a  natural  progression  of  increased 
complexity  (i.e.,  single  tooth  models  with  point  contact,  multi¬ 
tooth  models  with  point  contact,  multi-tooth  models  with 
deformable  contact,  etc.)  As  computers  became  more  powerful 
and  the  technology  became  more  mature,  eventually  entire 
planetaiy  gearsets  were  rolled  through  mesh  with  2D  models.  3D 
modeling  needs  to  go  through  a  similar  evolution. 

CONCLUSION 

Considerable  progress  has  been  made  in  automating  3D 
analysis  of  spiral  bevel  gears  in  mesh.  Programs  have  been 
written  that  start  with  machine  tool  settings  and  can  create  finite 
element  models  of  any  mesh  density.  This  is  done  by  numerically 
evaluating  the  kinematic  motion  of  the  manufacturing  process  for 
generated  spiral  bevel  gears  and  creating  input  data  for  a 
commercially  available  3D  solid  modeler. 

Contact  between  the  deformable  gear  teeth  is  modeled  with 
the  automatic  generation  of  non-penetration  constraints.  This  is 
done  by  contact  algorithm  in  a  commercially  available  finite 
element  code.  The  contact  algorithm  eliminates  the  problems  of 
mesh  distortion  associated  with  calculating  the  location  and 
orientation  of  gap  elements. 

Although  the  original  finite  element  analysis  was  done  on  a 
supercomputer  the  software  used  is  now  available  on  a  desktop 
PC. 
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Figure  2.  Coordinate  system  orientation  to  generate  a  right- 
hand  gear  surface  (c})^  =  0  shown  here) 
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Figure  6.  Non-penetrating  constraints  used  in  contact 
algorithm. 

Figure  4.  Distorted  gear  sur&ce  FE  mesh  required  for  gap 
element  orientation. 
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Figure  7.  Typical  contact  stress  contours  from  contact 
Figure  5.  Typical  contact  stress  contours  from  gap  element  algorithm  FE  model, 

model. 
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Figures  8a,  b,  c,  d,  e,  f.  Plot  of  contact  forces  for  pinion  segment  rolling  through  mesh 


Figure  9.  Typical  plot  of  contact  forces  on  two  tooth  gear  segment. 
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